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ABSTRACT 

We perform two-dimensional, Point-Spread-Function-convolved, bulge+disk decompositions in the 
g and r bandpasses on a sample of 1,123,718 galaxies from the Legacy area of the Sloan Digital Sky 
Survey Data Release Seven. Four different decomposition procedures are investigated which make 
improvements to sky background determinations and object deblending over the standard SDSS pro- 
cedures that lead to more robust structural parameters and integrated galaxy magnitudes and colors, 
especially in crowded environments. We use a set of science-based quality assurance metrics namely 
the disk luminosity-size relation, the galaxy color-magnitude diagram and the galaxy central (fiber) 
colors to show the robustness of our structural parameters. The best procedure utilizes simultaneous, 
two-bandpass decompositions. Bulge and disk photometric errors remain below 0.1 mag down to bulge 
and disk magnitudes of g ~ 19 and r ~ 18.5. We also use and compare three different galaxy fitting 
models: a pure Sersic model, a rtf, = 4 bulge -f- disk model and a Scrsic (free Ub) bulge + disk model. 
The most appropriate model for a given galaxy is determined by the _F-test probability. All three 
catalogs of measured structural parameters, rest-frame magnitudes and colors are publicly released 
here. These catalogs should provide an extensive comparison set for a wide range of observational 
and theoretical studies of galaxies. 
Subject headings: galaxies : fundamental parameters, galaxies : evolution 



1. INTRODUCTION 

The properties of disks and spheroids in the local 
Universe are the direct outcomes of the hierarchical 
assembly of galaxies over cosmological time. It is 
well known that the structural parameters of these im- 
portant galaxy sub-components obey scaling relations 
between size, lumin osity and internal velocity (e.g., 
iCourteau et al.l[2"007[ ). The amount of stellar mass found 
in disks and bulges places strong constraints on the 
galaxy me rger tree from ACDM N-body simulations 
(I Hopkins^ t al. 201( 3) and the secular forma tion of bulges 
(e.g.. lKormendv.,1979t IAthanassoulall200"3l) . Some prop- 
erties such as disk size can be traced to properties of host 
dark matter haloes such as specific angular momentum 
([Fall & Efstathiou 1980; Mo. Mao & White 1998), and 
more information is needed for a truly large sample of 
galaxies to improve the recipes used in semi-analytical 
approac hes to the treatment of t he so-called "sub-grid" 
physics (iGover nato et al.' '2004*: [Robertson et all 120041: 
iDutton fc vaiT dcn Bosch 2009). 

Thanks to the availability of large datasets and 
powerful computing clusters, it is now possible to 
analyze large samples of galaxies in a fully quantitative 
way. Quantitative measurements of galaxy properties 
have made s ignifica.nt progress oyer th e last few years 
jPeng et al.' '2002'; 'Simard et al." "2002"; 'de Souza ct al." 
2004^^ Lotz et al. 2004; Consclicc 2006; Pignatelli ct al. 
2006( 1 and now provide an important framework for 



comparisons between observation and theory. A number 
of previous works have looked at the morphologies 
of SDSS galaxi es using a wide range of subsamples 
and techniques ("Shcn e t al.l 12003 ': 'Blanton ct al.' ' 2005al: 
Kellv & McKav 2005; Be nson et a l. 2007; Fukuait a^etall 
2007; Lintott et al. 2008': iLa Barbera et al.l 120101 : 
Wijesinghc ct al. 2010). This paper presents the largest 
catalog of bulgc+disk structural parameters measured 
to date of galaxies in the Sl oan Digital Sky Survey Data 
Release Seven (SDSS DR7, lAbazajian eFaII[2009^ . This 
catalog has already been used to study visual versus 
quant itative morphologi e s (C h eng et al. 20 IQi), g alaxy 



pairs (jEUison et al.ll200l 120101: iPatton et al.ll201ll). disk 
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scaling r elations at low and hig h redshift (IDutton et afl 
l2011allbD . the disk size function (jKanwar et al.ll2008l and 
Simard 2011, in preparation), t he ev o lution of cluster 
galaxies with reds hift (Simard et al.l l2009l ). compact 
groups of galaxies (jMendel et al.ll201lD. and the disk 
and bulge stellar mass functions (jThaniavur fc Simard! 
1201 If ). The paper is organized as follows. The data 
and the bulge-|-disk decompositions are described in 
Sections [2] and |3l The catalog, science-based quality 
assessment metrics and a comparison between different 
fitting models are presented in Section (4] The cosmology 
adopted throughout this paper is (i?o,^m,^A) = (70 
km s-i Mpc-\ 0.3, 0.7). 

2. DATA 

The data come from the Legacy area of the SDSS. We 
selected a photometric sample from the PhotoPrimary 
table. The two main galaxy selection criteria were Pet- 
rosian magnitude (corrected for Galactic extinction) and 
morphological classification. We first selected objects 

with 14 < mpetro.r,corr < 18, whcrC mpetro,r.corr is the 

r-band Petrosian magnitude corrected for Galactic ex- 
tinction according to the extinction values given in the 
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SDSS database. We also selected extended objects with 
morphological type Type=3. In addition to these two 
criteria, we also required that the sum of the flags DE- 
BLENDED.AS.PSF and SATURATED be zero to ehmi- 
nate objects that were found to be unresolved children of 
their parents as well as saturated objects. The full query 
reads: 

select objid from drT.PhotoPrimary where 
flags & (dbo.fPhotoFlagsC SATURATED') + 
dbo.fPhotoFlags('DEBLENDEDJiS_PSF'))= and 
(petroMag_r-extinctionjr) between 14.0 sind 18.0 
and Type=3. 

This query returned 2,195,875 objects from the SEGUE 
and Legacy areas. In order to select only objects from 
the Legacy area, we used the DR7 sky coverage table 
allrunsdr7db . paiO to obtain the list of run/rerun com- 
binations in this area, and we then filtered the output of 
the query using the list. The total number of objects in 
our Legacy-area photometric sample is 1,123,718. 

Given the importance of the links between the 
morphological and spectroscopic properties of galax- 
ies, we also defined a spectroscopic galaxy sample by 
selecting objects from our photometric sample with 
'mpetro,r,corr ^ 17.77 and the spectrum of a galaxy as 
defined by the keyword SpecClass in the SpecPhoto 
database table. Specifically, we selected objects with 
SpecPhoto. SpecClass=2. These spectroscopic selection 
criteria yielded 674,701 galaxies. The faint limit is 
the completeness l imit of the SD SS Main Galaxy spec- 
tr oscopic sampl e (iStrauss et al.l 2002 1. As discussed 
in IStrauss et al.l (|2002D . the nominal surface brightness 
limit of the SDSS spectroscopic sample is r-band /iso.r = 
24.5 mag arcsec"^. However, objects in the range 23.0 
mag arcsec"^ < /iso.r £ 24.5 mag arcsec"^ were targeted 
only when the local and global sky values were within 
0.05 mag arcsec"^ of one another. We therefore set the 
faint surface brightness limit of our spectroscopic sam- 
ple to /X5o,r = 23.0 mag arcsec"^ to retain a complete 
sample. This limit further cuts the number of galaxies 
down to 671,425. The number of objects excluded by 
this cut is small as expected from Figure [T] Following 
iShen et aD (|2003f ). we then excluded a small number of 
galaxies with redshift z < 0.005 whose distances may be 
severely contaminated by their peculiar velocities. The 
final total number of objects satisfying our spectroscopic 
selection criteria was 670,131. 

g- and r-band images of all galaxies in our photometric 
sample were analyzed as described in Section [3l The red- 
shift distribution for the spectroscopic sample is shown in 
Figured] The bulk of the sample is at z ~ 0.1. The peaks 
in this distribution are due to real large-scale structures 
within the SDSS survey volume. 

3. ANALYSIS 

3.1. Preprocessing of SDSS Images 

Information for each galaxy to be analyzed was fetched 
directly from the SDSS catalog server using a remote 
data mining pipeline based on the Python program- 
ming language and the wget protocol. The fields 
objID, ra, dec, run, rerun, camcol, field, obj, 
petroMagjr, z, and extinction_r were extracted from 

* http://www.sdss.org/dr7/coverage/allrunsdr7db.par 
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Figure 1. Apparent r-band Petrosian magnitude and surface 
brightness distribution of galaxies in our photometric sample. Cor- 
rected for Galactic extinction using SDSS values. 
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Figure 2. Redshift distribution of galaxies in our spectroscopic 
sample (14 < mpetro,r,corr < 17.77, /^5o,r < 23 mag arcsec"^, 
SpecPhoto . SpecClass=2) . 

the PhotoPrimary table. The atlas (prefix "fpAt- 
las" ) , corrected (prefix "fpC" ) and point-spread-function 
(prefix "psField") images were then retrieved directly 
from the SDSS Data Archive Server (DAS). The loca- 
tions of these images for a given galaxy were given by 
run, rerun and csuncol. Galaxy positions (rowc and 
cole) on the corrected images were also extracted from 
PhotoPrimary. The SDSS corrected images are survey 
images that have been bias subtracted, fiat-fielded and 
purged of bright stars. 

Our bulge-f disk decomposition procedure requires in- 
put thumbnail images for each object. These thumb- 
nails were prepared from SDSS images using the fol- 
lowing pre-processing steps: (1) extraction of the atlas 
and point-spread-function images, (2) creation of the sci- 
ence thumbnail images, and (3) creation of the mask 
thumbnail image. We used the SDSS software utilities 
readAtlasImage43 and read_PSli3 for the atlas and PSF 

^ http://www.sdss.org/dr7/products/images/read_atlas.html 
^ http://www.sdss.org/dr7/products/images/read_psf.html 
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images respectively. readAtlasImages requires a band- 
pass and object ID as inputs while read_PSF can output 
a SDSS-provided PSF at any object position (given by 
rowc, cole) in a given bandpass. 

3.2. Defining Objects on SDSS images 

There is no unique procedure for deblending objects 
and defining their object-sky boundaries on astronomi- 
cal images. Two different approaches can be adopted to 
tackle this problem, and these two approaches are exem- 
plified by the SDSS PHOTO pipeline and the SExtractor 
software (jBertin fc Arnouta 1996). Both approaches are 
shown in Figure |31 In the case of SDSS deblending (Sec- 
tion [3?2TT]), the pipeline tries to isolate the flux in a given 
object to create a new image of what this object would 
look like if only this flux was shown. The object area 
of interest is then determined from this new deblended 
image. It is important to note that this object area may 
extend to areas that used to be occupied by neighbor- 
ing galaxies, and may affect the photometry if some flux 
from the neighboring galaxies that were removed is left 
behind. SExtractor deblending (Section [3. 2. 2p does not 
attempt to create a new image showing the deblended 
object. It uses a multi-threshold method to look at the 
flux tree of the object in its original image. As it moves 
up the tree going from one isophotal level to the next, it 
looks at the flux in each branch. If the flux in a branch 
is above a set minimum fraction of the total flux in the 
tree, then this branch is deemed to be a separate ob- 
ject altogether. Once all the objects have been iden- 
tified, SExtractor then produces a segmentation image 
that assigns pixels to the different objects. In the case 
of a galaxy with a close neighbor, the net result is a 
"ridge" of segmentation values at the saddle point be- 
tween the two adjacent objects. The segmentation of an 
object never extends over the area occupied by its neigh- 
bors (see upper left-hand panel of Figure |31). It is clear 
from Figure 13] that each object segmentation area misses 
a significant fraction of the flux in the object because 
it spills over into the segmentation area of its neighbor. 
The SExtractor magnitudes measured from these areas 
would therefore systematically underestimate the bright- 
nesses of the two objects. However, if we also flt some 
surface brightness model (e.g., bulge-f-disk) to the flux in- 
side a given object segmentation area and then integrate 
the flux in this model out to large radii, then the missing 
flux is recovered. As we will show later, such combina- 
tion of SExtractor deblending and model photometry can 
produce more reliable magnitudes and colors in crowded 
environments. 

3.2.1. SDSS Dehlendmg 

The atlas images provide object fluxes in pixels deemed 
to belo ng to the object by the SDSS deblending algo- 
rithm (jLupton et al.|[2001[ ). and they have zero object 
flux everywhere else. The SDSS pipeline adds a constant 
"softbias" level of 1000 counts to all pixels in the atlas 
and PSF images. This softbias level was subtracted from 
the images, so the atlas image sky pixels therefore have 
a value of zero. However, these atlas images cannot be 
used as is for our morphological analysis. The isophotal 
cut chosen to deflne an object is too high, and the wings 
of galaxy profiles where useful information can still be 



found have been truncated from these images as a result 
of this choice. Our science thumbnail images were there- 
fore produced by extending the object pixel area in the 
atlas images using pixels from the corrected image of the 
field in which the object resides. We ran SExtractor on 
the corrected image to identify sky pixels. Pixels from 
the corrected image were added to the science thumbnail 
only if they were identified as sky pixels. 

The mask image is a critical element of the analysis as 
it specifies the object pixel brightness distribution that 
will be "seen" by the fitting algorithm. Mask pixels could 
take on three different values: 1 (object), (sky) and 
—2 (excluded from fit). Initially, our main consideration 
was to remain consistent with the object deblending per- 
formed by the SDSS pipeline. We did not perform our 
own object deblending at first. We used the atlas image 
itself as a starting point for the creation of the mask im- 
age. For pixels in the mask image overlapping with the 
atlas image, we assigned them a value of 1 when the atlas 
image pixel value was greater than 1000 and a value of 
—2 otherwise. We then used the SExtractor segmenta- 
tion image of the corrected image to complete the mask. 
Pixels in the final mask were assigned a value of 1 if atlas 
mask pixel value was 1 and a value of if atlas mask was 
—2 and SExtractor mask value was 0. The net result was 
a mask with the same dimensions as the atlas image but 
with more sky pixels than the original atlas images that 
preserves the object deblending of the SDSS pipeline. 

3.2.2. SExtractor Dehlendmg 

The procedure for preparing GIM2D postage stamp 
images using SExtractor rather than SDSS deblend- 
ing is considerably simpler. It starts directly from the 
SDSS corrected image frame on which the object is lo- 
cated. SExtractor was run on this corrected frame to 
measure the parameters XJMAGE, YJMAGE, BACK- 
GROUND, and ISOAREA.IMAGE using a value of 
0.00005 for the SExtractor deblending contrast param- 
eter DEBLEND_MINCONT and a 1-cr isophote for the 
analysis area. The area extracted around each object for 
the GIM2D postage stamp images was set to five times 
the area given by ISOAREA.IMAGE. The GIM2D de- 
compositions were performed over all pixels fiagged as 
object or background in the SExtractor segmentation 
image. Objects in the segmentation images of the SDSS 
corrected frames are sharply delineated by the location 
of the isophote corresponding to the detection threshold 
because SExtractor considers all pixels below this thresh- 
old to be background pixels. However, precious informa- 
tion on the outer parts of the galaxy profile may be con- 
tained in the pixels below that threshold, and fits should 
therefore not be restricted only to object pixels to avoid 
throwing that information away (This is an analogous 
problem to the sharp isophotal limit in SDSS atlas im- 
ages). Pixels belonging to objects in the neighborhood of 
the primary object being fit are masked out of the fitting 
area using the SExtractor segmentation image. As noted 
earlier, the fiux from the primary object that would have 
been in those masked areas in the absence of neighbors 
is nonetheless properly included in the magnitude mea- 
surements because GIM2D magnitudes were obtained by 
integrating the best-fit models over all pixels. 

3.3. Sky Background Level Determination 
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Figure 3. The difference between SDSS and SExtractor deblending for object SDSS ObjID = 587729781199864370. This is the central 
galaxy in the image, and it has a close companion to the lower left. Top left: Object image from the SDSS corrected frame. Top right: 
Segmentation pixel mask produced by SExtractor as described in Section 13.2.21 The white region is the area defined as the object by the 
SExtractor deblending. Bottom left: Object i mage given by the SDSS deblender. Bottom right: Segmentation pixel mask constructed from 
the SDSS deblending as described in Section 13.2.11 The white region is the area defined as the object by the SDSS atlas image. Note how 
it now extends over the area where the close companion used to be. 



Systematic errors in the sky background level determi- 
nation dominate the errors on the measured structural 
parameters. An erroneous sky level can also mislead 
the bulge+disk decomposition algorithm into introduc- 
ing unphysical bulge or disk components. For example, 
a very large disk component would result from under- 
estimating the sky level because the positive sky offset 
would look like such a component to the algorithm. It is 
therefore critical to measure the best sky possible. 

We initially adopted the SDSS sky background levels 
for our bulge-|-disk decompositions for the sake of con- 
sistency. We used the sky levels given by the keyword 
SKY in the headers of the corrected images. The sky 
level for a given corrected frame was subtracted from all 
the GIM2D science postage stamp images of the objects 
extracted from this frame. The sky level was then fixed 
to zero for the bulge-f-disk decompositions. As discussed 
later, this procedure did not produce sky background 
levels that were good enough for our decompositions. 

We then used sky background levels and standard de- 
viations determined by GIM2D. GIM2D first uses ah the 
pixels in the science thumbnail image fiagged as back- 
ground pixels (flag value of zero) in the SExtractor seg- 
mentation image. GIM2D further prunes this sample 
of background pixels by excluding any background pixel 
that is closer than 4''0 from any (primary or neighbor- 
ing) object pixels. This buffer zone ensures that the flux 
from all SExtracted objects in the image below all the 
1.0-abkg isophotes does not significantly bias the mean 
background level upwards and artificially inflate abkg ■ A 
minimum of 20,000 sky pixels was imposed on the area 
of the sky region. In cases where the number of sky pix- 
els in the input science thumbnail image was insufficient, 
the original SDSS corrected image was searched for the 



20,000 sky pixels nearest to the object. Background pa- 
rameters were re-calculated with GIM2D before fitting, 
and the residual background levels were then frozen to 
their recalculated values for the bulge-|-disk fits. 

3.4. Two-Dimensional Bulge+Disk Decompositions 

Galaxy structural parameters were measured from 
bulge-t-disk decompositions per formed using versio n 3.2 
of the GIM2D software package (jSimard et al.ll2002[ ). We 
used the sum of a pure exponential disk and a de Vau- 
couleurs bulge (Sersic index n^ = 4) as our galaxy image 
model. The free fitting parameters of this model were the 
total fiux F in data units (DU), the bulge fraction B /T 
(= for pure disk systems), the bulge semi- major axis 
effective radius re, the bulge ellipticity e (e = 1 — 6/a, 
b = semi- minor axis, a = semi-major axis), the bulge po- 
sition angle of the major axis (pb on the image (clockwise, 
y-axis = 0), the semi- major axis exponential scale length 
rd (also denoted h in the literature), the disk inclination 
i (face-on = 0), the disk position angle 4)d on the image, 
and the dx and dy offsets of the model center with respect 
to the SDSS object position on the sky. The background 
residual level db and the bulge Sersic index Ub were held 
fixed for the fits. The position angles (pb and ipd were 
not forced to be equal for two reasons: (1) a large differ- 
ence between these position angles is often a signature of 
strongly barred galaxies, and (2) some observed galaxies 
do have bona fide bulges that are not quite aligned with 
the disk position angle. In addition to bulge-t-disk struc- 
tural parameters, GIM2D computes a variety of asym- 
metry indices that can be used to quantify galaxy sub- 
structure. The impact of substructures on GIM2D mea- 
surements has been extensively simulated and discussed 
in iS imard ct al., (|2002,) . We specifically chose to use a 
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bulge+disk model rather than a more complicated com- 
bination (e.g., bulge-l-disk-l-bar(-l-others)) for the sake of 
keeping our SDSS measurements consistent with what is 
achievable with galaxy structural measurements at high- 
redshift. It is already a challenge to perform reliable 
bulge+disk decompositions of high redshift galaxy im- 
ages even with the spatial resolution of the Hubble Space 
Telescope, and it is not currently possible to use struc- 
tural models with more components without compromis- 
ing convergence and avoiding parameter degeneracies. 

Fits in the g- and r- bands were done using the two- 
bandpass se parate and si r nultan eous fitting procedures 
described in lSimard et al.l ()2002D . For the separate, two- 
bandpass fits, the g- and r-band images were fitted inde- 
pendently. For the simultaneous, two-bandpass fits, the 
images were fitted by forcing the bulge radius, ellipticity 
and position angle and disk scale length, inclination and 
position angle to take on the same values in both band- 
passes. The main advantage of the simultaneous fits is to 
minimize the errors on the measured structural parame- 
ters because all of the data is used at once. Their main 
disadvantage is that they are blind to color gradients, 
and significant color gradients would make this approach 
invalid because sizes would be significantly different in 
different bandpasses. To test the validity of the applica- 
tion of simultaneous fitting to the SDSS data, we plotted 
the ratio of the g-band and r-band disk scale lengths from 
separate fits as a function of disk magnitude (Figure |4]). 
The two disk scale lengths are within a few percent of 
one another and in agreement wit h a recent cosmolog i- 
cally based disk formation model (jDutton et al.ll20l"lal ). 
As also shown in Figure 21 the bulge effective radii mea- 
sured separately in the two bandpasses agree to 2%. We 
therefore deem the simultaneous fitting approach to be 
valid. There is one important difference betweenjthe si- 
multaneous fitting p rocedure of|S imard ct al. (2002!) and 
the one used here. iSimard et al.. (,2002.) used the same 
dx and dy offsets in both bands because their images 
had been registered to the same pixel grid in both bands. 
This is not the case with the SDSS images, and using the 
same set of offsets produced "positive-negative" residu- 
als in the cores of many galaxies. We therefore let dx 
and dy be different in g and r, i.e, we let {dx)g, {dy)r, 
{dx)r and {dy)r be free parameters in the simultaneous 
fits. The simultaneous fits therefore used a total number 
of fourteen free parameters {Fg, Fr, {B/T)g, {B/T)r, r^, 
e, 0b, rd, i, (t>d, {dx)g, {dx)r, {dy)g and {dy)r). 

Four sets of GIM2D fits were performed on the 
spectroscopic sample defined in Section [5] with dif- 
ferent combinations of object deblending, sky level 
determinations and fitting procedures: (1) Sepa- 
rate fitting + SDSS object deblending + SDSS 
sky level ("SEP+SDSSDEBL+SDSSBKG"), (2) Sep- 
arate fitting + SDSS deblending -I- GIM2D sky 
level ("SEP-FSDSSDEBL+GM2DSKY"), (3) Simulta- 
neous fitting + SDSS deblending -I- GIM2D sky level 
("SIM-fSDSSDEBL+GM2DSKY"), and (4) Simultane- 
ous fitting + SExtractor deblending + GIM2D sky level 
("SIM+SEXTDEBL+GM2DSKY"). This last set of fits 
was adopted as our preferred set on the basis of quality 
assessment tests described in Section HTTl and the proce- 
dure was then run on our full photometric sample of 1.12 
million galaxies. Out of all the galaxies in this sample, 
required files could not be obtained from the SDSS im- 



age server for 58 galaxies, and GIM2D convergence was 
not achieved for 7361 objects (failure rate of 0.66%). Vi- 
sual inspection of these failures traced most of them back 
to artefacts (e.g., image defects, bright star spikes, bits 
of a "shredded" galaxy) in the SDSS catalog. Figure [S] 
shows a detailed example of a bulge-l-disk decomposition 
of a single galaxy, and the summary decompositions for 
twenty galaxies drawn at random from our sample are 
shown in Figure [6] 

In addition to our canonical ri;, = 4 bulge -t- disk fit- 
ting model (denoted as "n4" hereafter), we used two 
other fitting models for comparison. The first additional 
model was a free Ub bulge -t- disk model (denoted as 
"fn" hereafter) with rib allowed to vary from 0.5 to 8, 
and the second one was a single component, pure Sersic 
model (denoted as "pS" hereafter). The fitting param- 
eters of the latter were Fg, Fr, r^, e, 4>b, {dx)g, {dx)r, 
{dy)g, {dy)r and Ug where Ug is now the global galaxy 
Sersic index, and Ug was also allowed to vary from 0.5 
to 8. All of these additional fits were performed with 
the same object deblending and sky background level as 
the "SIM-l-SEXTDEBL-f GM2DSKY" dataset described 
above. Structural parameters from all three fitting mod- 
els are compared in Section 14.21 and given in the data 
tables in Section H^ 

All of our galaxy model fluxes were converted to SDSS 
magnitudes using: 



rrigaiaxy = -2.51ogio(F/i) - xsecz - uiq 



(1) 



where F is the model flux, t is the exposure time of the 
SDSS corrected images (53.907456 seconds), sec z is the 
airmass, x is the extinction coefficient, and ttt-o is the 
magnitude zeropoint. Values for these coefficients were 
taken directly from the SDSS database table Field and 
applied on a field- by-field basis. The photometric coeffi- 
cients from this table put our fiux measurements on the 
SDSS ugriz magnitude system. We did not transform 
our magnitudes from the SDSS to the AB system. The 
largest offset between the two systems is found in the 
u— band where uab = usdss — 0.04. The offsets in the 
other bands are all 0.01 mag or less. 



3.5. Rest-Frame Quantities 

All angular sizes were converted to physical sizes ac- 
cording to the equation: 



R^r- 



dz' 



1000Hol + 2 7o v/"™(l + z'f + Qa 



(2) 



where z is the redshift, c is the speed of light, r is the 
measured angular size in radians, and R is the corre- 
sponding physical size in kiloparsecs. Equation [5] is only 
valid for flat {ilk = 1 ^ f^m — ^A — 0) cosmologies. 

The g-band and r-band galaxy, bulge and disk rest- 
frame absolute magnitudes were computed from the 
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Figure 4. Left-hand panel: Log ratio of g- to r-band disk scale lengths measured from separate bandpass fits as a function of disk 
magnitude for galaxies with r-band bulge fraction (B/T)r < 0.5 and apparent half-light radius > 2"0. The average ratio of the two disk 
scale lengths is lO''"^ = 1.05. Right-hand panel: Log ratio of g- to r-band bulge effective radii measured from separate bandpass fits as a 
function of bulge magnitude for galaxies with r-band bulge fraction {B/T)r > 0.5 and apparent half-light radius > 2''0. The average ratio 
of the two bulge effective radii is lO"'"^ = 1.02. 



GIM2D apparent magnitudes according to the equations: 
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where mg^g and rrir^g are the y-band and r-band observed 
GIM2D galaxy model magnitudes from Equation [I] eg 
and Cj. are the line-of-sight galactic extinctions in mag- 
nitude from the SDSS database, {B/T)g and {B/T)r are 
the GIM2D observed bulge fractions, DM{z) is the dis- 
tance modulus for redshift z, and kg and kr are the 
fc-corrections. The fc -corrections were comput ed using 
kcorrect version 4 (jBlanton fc Roweii l2007f) . Rest- 
frame magnitude s ar e included in the data table pre- 
sented in Section l4?3l 



3.6. Sample Volume Corrections 

Volume corrections must be applied to galaxies in our 
sample to properly account for the effects of selection on 
the visibility of different classes of galaxies. The SDSS 
spectroscopic sample is complete down to a r-band mag- 
nitude of 17.77 and an effective surface brightness ^50, r of 
23 mag arcsec"^ (corrected for Galactic extin ction). Our 
Vmnx- corrections are similar to the ones in iShen et alJ 
(j2003| ). First, the magnitude range rmin ^ r < r^ax 
corresponds to a maximum redshift Zmax,m and a mini- 
mum redshift Zmin,m' 

dL{Zmax,m)^dL{z)10-"-^^''-^'--^^ (4a) 

rfL(z„,„,™) = (iL(z)10-"-2('^-— ) (4b) 



where d^^z) is the luminosity distance at redshift z. 
The surface brightness limit constrains Vmax for a given 
galaxy mainly through the (1 -\- z)'^ dimming effect. The 
maximum redshift at which a galaxy of surface brightness 
M50,r at z can still be detected within a limiting surface 
brightness ^um = 23 is given by: 



^max.fi 



= (l + z)10(23-0-''--)/l"-l 



(5) 



One difference between the Vm 



corrections in 



iShen et all (|2003| ) and the ones used here is that we do 
not impose a minimum size limit to galaxies in our sam- 
ple. The real maximum and minimum redshifts, z^ax 
and Zmim for a given galaxy are therefore: 



: max(zi 
: min(z; 



vmn.ra 1 U.UUO j 
Zmnx-^J 



max.ini '^max^. 
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With the above limits, Vmax is then given by the inte- 
gral: 



K, 



47r 



dnf{0,(i)) 



.(9,0) 



d\{z) 



H{z){l + z) 



cdz{7) 



where H{z) is the Hubble parameter at redshift z, c is 
the speed of light, f{9, (f) is the sampling fraction as a 
function of position on the sky, and fi is the solid an- 
gle. We take the function f{9, </>) to be a constant over 
the entire 8032 square degree areal coverage of the SDSS 
DR7 Legacy survey. The angular part of Equation[7]thus 
becomes (l/47r) (8032/41253) (47r) = 0.1947. V^ax values 
are given in the data table presented in Section 14.31 



4. RESULTS 

4.1. Quality Assessment 

We based the quality assessment of our measured 
structural parameters on three science-motivated met- 
rics. 
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Figure 5. Example of a GIM2D bulge+disk decomposition (SDSS objID = 587741708861964592) . From left to right: Galaxy cutout 
from SDSS corrected image, GIM2D output model image and GIM2D residual image. The top two rows show the bulge component only, 
the middle two rows show the disk component only, and the bottom two rows show the full bulge+disk model. The labels at the bottom 
of the bulge component images are effective radius in kiloparsecs, elliptieity (= for circular) and position angle (measured clockwise from 
+y axis. The labels at the bottom of the disk component images are disk scale length in kiloparsecs, inclination angle (= for face-on) 
and position angle (measured clockwise from +y axis). The labels are the bottom of the bulge+disk model images are rest-frame absolute 
magnitude, radius in kiloparsecs and bulge fraction. The minimum and maximum values for the greyscale stretch are <bkg> — 5o'bkg ^^nd 
<bkg> +30cr)3kg respectively where <bkg> and (Tbkg ^^re the mean and dispersion of the background pixel values. 



4.1.1. Size-Lummosity Relation of Disks 

Figure [7] shows the size-luminosity of disks in the g- 
band. Disks are expected to follow a well-defined and 
well-known scaling relation, and most of the disks do 
follow such a relation. However, some disks in the 
"SEP+SDSSDEBL+SDSSBKG" photometric dataset lie 
on an "upper branch" of large disks that appear to be 
2-3 times larger at a given disk luminosity. Plotting 
the same relation for the M-band showed an increasing 
bimodality in the sizes of disks. Visual inspection of 
these apparently large disks showed that they were in 



fact caused by constant sky background residuals left 
after subtracting the SDSS sky values. As explained 
in Section 13.31 our initial procedure assumed that the 
SDSS sky values were "perfect" and that the GIM2D 
fits could be performed by fixing the sky value to zero 
once the SDSS sky values had been subtracted from the 
input science images. These constant positive sky resid- 
uals were being fit by the algorithm with a very large 
(and thus very flat) disk component. The false disks 
lie on a linear relation because of the strong covariance 
between luminosity and size errors through the relation 
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Figure 6. Mosaic of r-band GIM2D bulge+disk decompositions. Three postage stamp images are shown for each galaxy: galaxy cutout 
from SDSS corrected image, galaxy bulge+disk GIM2D model convolved with its point-spread-function, and the GIM2D residual image. 
Each corrected image cutout is labelled by objID, apparent GIM2D model magnitude, GIM2D model half-light radius in arcsecs, and its 
SDSS redshift. Each GIM2D model image is labelled by GIM2D model rest-frame absolute magnitude, GIM2D model half-light radius in 
kiloparcsecs and GIM2D model bul ge fraction. E ach GIM2D residual image is labelled by three asymmetry indices: _RT1_2, RA1_2 and 
S2 = RT1J2 + RA1_2 (as defined in I Simard et al . (2002, 2009)). The minimum and maximum values for the greyscale stretch for a given 
galaxy are <bkg> — 5crbkg and <bkg> +30crbkg respectively where <bkg> and crbkg are the mean and dispersion of the background pixel 
values. 
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Figure 7. Top panel: g-band disk size-luminosity rela- 

tion without volume corrections for face-on (0.75 < b/a < 
1.0) disks from separate fits with SDSS background levels 
("SEP-l-SDSSDEBL-fSDSSBKG"). Bottom panel: Same disk 
sample but with GIM2D-determined sky background levels 
("SEP-I-SDSSDEBL-I-GM2DBKG"). Note the disappearance of 
the "large disks" branch with the use of GIM2D sky levels. 



F — 27rEo'"d for the total flux F of an exponential disk 
with central surface brightness So and scale length rd- 
Recomputing local sky background levels following the 
procedure described in Section 13.31 corrected this prob- 
lem as shown by the "SEP-hSDSSDEBL-f GM2DBKG" 
photometric dataset in the bottom panel of Figure [T] 
The u-band relation was also corrected by the same pro- 
cedure. The disk size-luminosity relation is therefore a 
sensitive test of our sky background level determinations. 
It is worth noting here that this sky background prob- 
lem is different from the well-known p roblem with the 
photometry of bright galaxies in SDSS (" Abazaiian et al.l 
[2009. and references therein). The SDSS photometric 
pipeline produces systematic errors in the estimation of 
the sky near bright (r < 16) galaxies that cause their 
fluxes and scale sizes to be underestimated. The sky 
background problem encountered here was related to 
positive sky residuals. These effects of these residuals be- 
came especially significant here because sky level errors 
domi nate the systematic errors in bulge-f disk decomposi- 
tions (jSimard et al.l[2003 ) . It will also be more significant 
in more crowded environments (including close pairs of 
galaxies). 
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4.1.2. Integrated Galaxy Magnitudes and Colors 

Figure [8] compares our GIM2D-based total magni- 
tudes with the SDSS Petrosian and model magnitudes 
as a function of bulge fraction. There are two differ- 
ent types of SDSS model magnitudes. The first type 
(called "modelMag" ) comes from fitting the object with 
a pure deVaucouleurs and a pure exponential profile, 
determining which profile fits the object better and 
computing the magnitude of this profile. The "model- 
Mag" magnitudes are widely used throughout the liter- 
ature because they are the ones explicitly given in the 
SDSS database. The second type of model magnitudes 
(called "cmodel") is computed from the composite flux 
given by the linear combination of the two profiles, i.e., 
Fcompostte = {fracDeV)FdeV + (1 - fracDeV)Fexp (the 
value of fracDev is stored in the SDSS database). Al- 
though it would be more fair to compare the GIM2D 
magnitudes to the "cmodel" magnitudes, we decided to 
compare our photometry with the most widely used set 
of SDSS model magnitudes to emphasize that the use of 
SDSS model magnitudes requires some caution for galax- 
ies of intermediate morphological types. This is obvious 
from looking at the panels in the first column of Figure |8l 
These panels do not depend on any GIM2D photometry. 
They show a "bifurcation" at intermediate bulge frac- 
tion values because the "modelMag" model magnitudes 
can only come from either a pure disk or a pure bulge 
model, and some objects with a mixture of bulge and 
disk components become ambiguous enough to such a 
binary classification that they end up in different bins. 

The second and third columns of Figure [H] show that 
the offset between GIM2D and SDSS magnitudes is in- 
dependent of bulge fraction. This lack of dependence 
on the type of galaxy light profiles mak es it easier to 
compa re the different sets of magnitudes. iBlanton et al.l 
(|2003f ) reported a systematic trend between their Sersic- 
based magnitude and SDSS Petrosian magnitudes as a 
function of the Sersic index value with Petrosian magni- 
tudes underestimating the g- and r-band fluxes by 20% 
in their high (n ~ 4) Sersic index galaxies. The bifur- 
cation at intermediate bulge fractions is also visible in 
the comparison between GIM2D and SDSS model mag- 
nitudes. 

Object deblending can have a significant impact on 
galaxy photometry especially in crowded environments 
(Section 13. 2p . The simplest case of a crowded environ- 
ment is a close pair of galaxies, and we teste d our GIM2D 
photo metry using the galaxy pair sample of lPatton et al.l 
(|2011f ). This galaxy pair sample also comes with a con- 
trol sample matched in stellar mass and redshift, and 
we used this control sample for a comparison between 
crowded and non-crowded environments. The control 
sample provides a cleaner comparison than the entire 
SDSS sample in the sense that selection effects due to 
galaxy properties such as luminosity and morphological 
type and technical considerations such as fiber collision 
are the same for both samples. Figures [9] and \W\ show 
color-magnitude diagrams for the Patton et al. galaxy 
pair and control samples. These diagrams are shown 
for different types of photometry. The GIM2D separate 
fits ("SEP+SDSSDEBL+SDSSBKG") produce a red se- 
quence and a blue cloud that are less well-defined than 
those based on Petrosian and model magnitudes. GIM2D 
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magnitudes are total magnitudes in the sense that the 
fluxes in the best-fit bulge and disk components are in- 
tegrated to infinity, and as such, they will have a lower 
signal-to-noise ratio than Petrosian magnitudes, say, that 
aim to minimize the amount of missing flux while opti- 
mizing the signal-to-noise ratio of the photometric mea- 
surements. The GIM2D fits also yield a large num- 
ber of bright galaxies over the range — 24 < Mr < 21. 
As discussed in Section 14.1.11 the GIM2D-determined 
sky backgrounds yielded better disk sizes and luminosi- 
ties. Adopting this new background level determination 
method also decreases the number of those bright out- 
liers, but significant scatter towards bluer and redder col- 
ors remains. 

Switching from bulge-t-disk fits performed separately in 
each filter to the ones performed simultaneously in both 
bandpasses immediately improves the integrated magni- 
tudes in two important ways. The number of blue and 
red outliers is now significantly lower, and the red se- 
quence and blue cloud are now as well-defined as for 
the SDSS photometry. The simultaneous fits produce 
smaller errors because they simultaneously make use of 
all the data available for a given object. Visual inspection 
of the remaining outliers showed them to be in crowded 
environments where the SDSS object deblending was fail- 
ing. Adopting the SExtractor deblending (in addition to 
the simultaneous fits) described in Section 13.2.21 elimi- 
nated nearly all of the remaining outliers to produce a 
very clean color-magnitude diagram. The lack of red out- 
liers in the final GIM2D color-magn itude diagrani i s espe - 
cially noteworthy. As explained in iPatton et al.l ()2011[ ) , 
the presence of these outliers in the official SDSS pho- 
tometry was ascribed in the literature to a new popu- 
lation of e xtremely-red pair galaxies (e.g., lAlonso et al.l 
I2OOI IPere z ct al. 2009; Darg et al. 2010). The GIM2D 
photometry given here shows that they are in fact due to 
crowding errors from incorrect deblending. The crowd- 
ing errors in SDSS preferentially scatter galaxy colors 
towards the red for the following reason. The r-band im- 
ages go deeper in surface brightness than the g-hand ones 
because the SDSS detectors are more red sensitive. The 
outer isophotes of two neighboring galaxies will therefore 
overlap more in the red than in the blue, and the SDSS 
deblender will leave more faint isophotal light behind in 
the red. When this residual light for a galaxy is included 
in the photometry of its neighbor(s), it yields (a) redder 
neighbor galaxy color (s). 

4.1.3. Fiber Colors 

Our third data quality metric is based on fiber appar- 
ent magnitudes and colors. The fiber magnitude of a 
given object is computed from the flux measured within 
the aperture of a spectroscopic flber. The diameter of 
this aperture is 3". The fiber magnitudes are given in 
the SDSS database. Given the apparent sizes of the 
SDSS galaxies with respect to this fiber aperture, fiber 
magnitudes (and colors) are most often measurements 
of the central part of a galaxy. The typical fiber cover- 
ing fraction, which is the ratio of the q-ban d Petrosian 
to fiber fluxes, is 30% (jEllison et al.ll200l . We can 
also compute fiber magnitudes from the GIM2D out- 
put model images for comparison with the SDSS fiber 
quantities. The fiber aperture on the GIM2D model im- 
ages was centered on the SDSS position of the galaxy 



on the sky {{dx,dy) — (0,0)) and not on the center 
of the GIM2D model itself which was allowed to move 
for the fits. Using the SDSS and GIM2D fiber mag- 
nitudes, we then computed a A(fiber color) given by 
{g-r)gim2d., fiber " {g - r) SDSS, fiber ■ This A(fiber color) 
is plotted in Figure [TT] as a function of projected galaxy- 
galaxy se2aration_on the sk y for galaxies in the pair sam- 
ple of IPatton erail ()201lD . The median A (fiber color) 
value is not zero due to PSF effects. The SDSS fiber 
magnitudes were computed by uniformly convolving all 
the SDSS images to a 2" seeing whereas the typical size of 
the SDSS PSFs used for the GIM2D model images is 1 '.'4 
FWHM. This narrower PSF FWHM spreads less bulge 
(i.e., red) hght outside of the fiber aperture thus making 
the GIM2D fiber colors slightly redder. A large positive 
or negative A (fiber color) indicates that a bulge-t-disk 
decomposition failed because it did not even reproduce 
the central parts of the galaxy. SDSS deblending pro- 
duces a large number of outliers, and the number of these 
outliers increases with decreasing pair separation. The 
number of galaxies with large A (fiber color) values in- 
creases significantly at separations less than 20 kilopar- 
secs which corresponds to angular separations less than 
10 "8 at z = 0.1. Switching to SExtractor deblending 
eliminates the vast majority of the outliers, and no trend 
is seen anymore as a function of pair separation. This is 
another example where SExtractor deblending produces 
more reliable photometry in crowded environments. 

4.2. Galaxy Fitting Models 
4.2.1. Selecting the Appropriate Model with the F-Statistic 

It is not always true that a galaxy will be appropriately 
modelled by a multi-component photometric model; in 
instances where galaxies are best represented by a pure 
bulge or disk or, alternatively, the data cannot support a 
more complex model decomposition due to low signal-to- 
noise or poor spatial sampling, it is unclear how to treat 
the output of our "n4" or "fn" fits. We therefore want a 
way to assess the appropriateness of different model de- 
compositions, and a means of deciding when a complex 
model, e.g. bulge-l-disk, is to be preferred over a pure 
Sersic profile. We approach this question of finding the 
appropriate fitting model using the i^-test to compare 
the x^ re siduals of our "pS" , "n4" , and "fn" fits (see Sec- 
tion [33]). The F-test provides a quantitative means of 
judging the relative likelihood of different structural de- 
compositions. Th ere are of course other approaches (e.g., 
I Allen et all I2006D . but the F-test seems to be a pow- 
erful and straightforward approach as we show below. 
For each model we estimate the reduced x^ using the fit 
residuals and the number of resolution elements per fit. 



n■ccs^ in this instance estimated as n^ 



"-pixGls/(7I'6^) 



where npixeis is the number of unmasked (object) pixels 
used in the fit and is the r-band seeing half-width at 
half maximum (HWHM) in pixels. The computed F val- 
ues for a given pair of models, e.g. "fn" vs. "pS", can 
then be translated into a corresponding probability that 
a smaller F would be observed under the null hypothesis 
that a galaxy is similarly well fit by the two decomposi- 
tions (i.e. that their rms flux residuals are comparable); 
these probabilities are given in Tables [1] and [21 In a 
given table, we only provide the probability relative to 
less complex models, such that the probabilities always 
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Figure 8. Comparison of r-band GIM2D magnitudes with Petrosian and model SDSS magnitudes ("modelMag") as a function 
of r-band bulge fraction. First column: SDSS Petrosian versus SDSS model. Second column: SDSS Petrosian versus GIM2D 
("SIM-I-SEXTDBL+GM2DBKG"). Third column: SDSS model versus GIM2D ("SIM+SEXTDBL+GM2DBKG"). Fourth column: SDSS 
model versus GIM2D ("SEP+SDSSDBL+SDSSBKG" denoted 'mgim2d,aep here). Note that the panels in the first column do not depend 
on any GIM2D photometry. 



correspond to the lower-tail probability of the appropri- 
ate i^-distribution. In other words, Pps is the probability 
that a bulge+disk model is not required (compared to a 
pure Sersic model) , and P„4 is the probability that a free 
nb bulge -|- disk model is not required (compared to a 
fixed ni, — 4 bulge -I- disk model). 

The distributions of F-test probabilities versus bulge 
fraction are shown in Figure I12i and they bring forth 
the limitations of the SDSS imaging in terms of spatial 
resolution and signal-to-noise. A bulge-|-disk model is 
clearly required for galaxies with 0.2 < {B/T)fn < 0.45, 
and galaxies with {B/T)fn > 0.75 do not require a 
bulge-(-disk model to fit their light profiles. If we set 
a value of Pps < 0.32 as the (Icr) threshold below which 
galaxies are likely to be genuine bulge-|-disk systems, 
then the fraction of these galaxies in our entire SDSS 
sample is 26%. The quality of the SDSS imaging is in- 
sufficient to determine bulge Sersic indices for galaxies in 
our selected range of apparent magnitudes as P„4 versus 



bulge fraction does not show any statistically significant 
differences between rif, = 4 and free Ub models. Only 9% 
of the galaxies have P„4 < 0.32. 

Figure [13] illustrates the usefulness of the F-test prob- 
abilities to select genuine bulge-|-disk systems. The axial 
ratio distribution of a sample of disks randomly inclined 
in space should be uniform between zero and one modulo 
some perturbations due to dust and/or bars. Looking 
at the top row of Figure [T^l one can see that the ob- 
served distribution is indeed uniform at low bulge frac- 
tion, but that it also tends towards the same axial ratio 
distribution as bulges towards higher bulge fractions if 
no -F-test selection is made. This behavior would argue 
that "disks" in highly bulge-dominated galaxies are not 
in fact real but were rather introduced by the fitting al- 
gorithm as an additional degree of freedom to model the 
outer wings of a single component galaxy. If we select 
only galaxies with Pps < 0.32, then one can see that the 
resulting disk axial ratio remains uniform even at the 
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Figure 9. Comparison of rest-frame color-magnitude diagrams for thc'Patton et al.' C2011') DR7 galaxy pair sample of (7V=22565): Top 
left: SDSS Petrosian magnitudes Top right: SDSS model magnitudes Middle left: GIM2D separate fits with SDSS sky background and 
deblending, Middle right: GIM2D separate fits with GIM2D sky background and SDSS deblending Bottom left: GIM2D simultaneous fits 
with GIM2D sky background and SDSS deblending. Bottom right: GIM2D simultaneous fits with GIM2D sky background and SExtractor 
deblending. The same greyscale was used in all subpanels. The bin size for the greyscale was {A(g — r), AMr) = (0.1, 0.05). The greyscale 
was replaced by the actual data points wherever the number of points in a given bin was less than 2. 
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Figure 10. Comparison of rest-frame color-magnitude diagrams for the'Patton et al." f20ir) DR7 galaxy pair control sample (N=290090): 
(Top left) SDSS Petrosian magnitudes, (top right) SDSS model magnitudes, (middle left) GIM2D separate fits with SDSS sky background 
and deblending, (middle right) GIM2D separate fits with GIM2D sky background and SDSS deblending, (bottom left) GIM2D simultaneous 
fits with GIM2D sky background and SDSS deblending, and (bottom right) GIM2D simultaneous fits with GIM2D sky background and 
SExtractor deblending. The same greyscale was used in all subpanels. The bin size for the greyscale was (A(g — r), AMr) = (0.1, 0.05). 
The greyscale was replaced by the actual data points wherever the number of points in a given bin was less than 2. 
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Figure 11. A(fiber color) as a function of galaxy pair separation in kiloparsecs with SDSS deblending {left-hand panel) and SExtractor 
deblcnding {right-hand panel). A{fiber color) is defined as {g — r)gi^2d, fiber ~ {d ~ ^)sDSS, fiber- The lines show the median fiber color 
(solid line, median value = 0.0182) from SDSS deblending and the 1-cr envelope (dashed lines, cr = 0,116), The same lines are reproduced 
in both panels. Note the excess of outliers in the left panel at small separations due to poor deblending. 



highest bulge fractions. We are thus able to select real 
disks in highly bulge dominated galaxies. 

4.2.2. Structural Parameter Comparison 

A detailed comparison between our three different fit- 
ting models could easily stand on its own as a separate 
paper, so we briefly comment here on the salient points 
of this comparison. We start by looking at measure- 
ments of the global galaxy structure from the different 
models namely the galaxy Sersic index ng and the bulge 
fraction B/T. The distribution of galaxy Sersic index 
values in Figure [14] shows a peak a sharp peak at Ug = 
0.5, a broader and larger peak at n^ = 1, no peak at rig 
— 4 and a peak at the maximum allowed Sersic value 
of 8. We will return to the peaks at rig = 0.5 and rig 
= 8 later. The large peak at ng = 1 reflects the fact 
that the local galaxy population is dominated by disk 
galaxies. The lack of a strong peak at rig — A is quite 
interesting because it argues for a lack of global struc- 
tural similarity in galaxies that are not disk-dominated. 
In terms of bulge fraction, measured values from both 
the rib — A and free n^ fits (denoted as {B/T)n4 and 
{B/T)fn respectively here) are well correlated with one 
another and with rig (Figure [M]) . {B/T)n4 and (B/r)/„ 
increase with rig and reach a value of one at n^ = 4. 
Two "branches" are seen beyond rig = 4: one at B/T ~ 
1 and another one at B/T ~ 0.5. The first branch is 
not surprising, but the second one may point to a poten- 
tial single- versus double-component degeneracy in the 
bulge-|-disk decompositions. In order to understand the 
behavior of {B/T) at high rig, we looked at the bulge and 
disk ellipticity difference Ae = e^ — e^ versus the ratio 
of the bulge and disk half-light radii Rb/d = ''e/Cl-SIrd) 
in the range 7 < n^ < 8 for two ranges in {B/T)n4- 
0.4 < (B/T)„4 < 0.7 and 0.9 < {B/T)n4 < 1.0. Ae 
is zero which means that the bulges and disks have the 
same ellipticities i.e., these galaxies may actually have a 
single component but the fitting algorithm may be using 
a bulge and disk components to model something that 
it cannot do even with (B/T)n4 = 1 for rib — 4. The 
dichotomy actually comes from the ratio of the radii. 
Galaxies with 0.4 < {B/T)nA < 0.7 all have Rb/d values 
of 0.13 with a dispersion of 0.05 whereas galaxies with 
0.9 < {B/T)n4 < 1.0 have values around 0.56 with a dis- 



persion of 0.16. There are six times more galaxies with 
0.4 < {B/T)n4 < 0.7 than with 0.9 < {B/T)n4 < 1.0. 
The peak in Rb/d where a galaxy ends up seems to de- 
pend on the spread in Ae. The dispersion in Ae for 
galaxies with 0.4 < {B/T)nA < 0.7 is twice the disper- 
sion in Ae for galaxies with 0.9 < {B/T)n4 < 1.0. A 
smaller Ae value for a given galaxy makes it more likely 
that the fitting algorithm will tend towards a single com- 
ponent model rather than a two-component model be- 
cause the algorithm will need two components to model 
a change of ellipticity with radius that does not come 
from PSF smearing. There is no significant difference 
between {B/T)n4 and {B/T)fn as a function of bulge 
Sersic index except at very low values of rib — 0.5 where 
{B/T)n4 can be considerably lower than {B/T)fn (Fig- 
ure I14p . This behavior at low Ub is expected given that 
trying to fit a bulge with a very flat proflle using a rel- 
atively peaky nb=4 component will force the algorithm 
to minimize the contribution of this nb=4 component 
as much as possible by converging t o a very low bulge 
fracti on. Some previous studies f e.g.. I Graham fc WorlevI 
|2008[ ) have reported a dependence of {B/T) on rib- Such 
a dependence is not seen here, but this may be due to the 
lack of constraint on bulge profile shape from the SDSS 
images of the galaxies in our sample as discussed next. 

The choice of a Sersic index value for the galaxy 
bulge profile has long b een deba ted in the literature 
(lAndredakis et all Il995t Ide Jong I [1.9 96; Bal cells et al 
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20101 : [Laurik ainen et al. 201Q), and our very large sample 
of free rib + disk decompositions should at first glance 
offer some insight on this interesting question. How- 
ever, the vast majority of the galaxies in our sample 
have images that do not have the required spatial resolu- 
tion and/or signal-to-noise ratio (Section 14.2.11 and Fig- 
ure [12]). If we select galaxies for which both Pps and P„4 
are less or equal to 0.32, then we obtain a subsample of al- 
most 53,000 galaxies for which we have good enough im- 
ages to study their bulge profile shape. Figure [T51 shows 
the distribution of n^, for this subsample. There are three 
important features to examine here. First, there is a peak 
at Ub — 0.5. We examined the distribution of galax- 
ies with 0.5 < rib < 0.55 in apparent bulge size Tg and 
bulge ellipticity e, and this distribution showed that es- 
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Figure 13. Bulge and disk axial ratio distributions for ri;, = 4 -|- disk models. In each panel, shaded and open histograms show the 
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b/a > 0.5 disk components at high B/T. Bottom row. Axial ratio distribution for the subsample of n;, = 4 -|- disk models where PpS < 0.32, 
i.e., where galaxies are likely to be genuine bulgc+disk systems. Note the significant reduction of low-inclination disks, particularly at 
B/T > 0.5. 



sentially all of these galaxies (TV ^ 6,700) were located 
in a peak at e ~ 0.7 and r^ ^ 1" -2" , i.e., 3-5 pixels. 
A combination of low nt and high e values is expected 
when the fitting algorithm tries to make the bulge pro- 
file as flat and as elongated as possible to try to fit a 
bar or include off-center components (point sources, very 
close mergers, etc.). We visually inspected a subsample 
of galaxies in this peak using the SDSS SkyServer (Ob- 
ject ExploreiQ to confirm this expectation. Second, the 
H}, distribution has a broad bump around ni, — 5.5 — 6.0. 
The de Vaucouleurs value of Tifc = 4 is not preferred for 
this subsample. However, choosing n{, = 4 for the entire 
sample is still a reasonable choice for the following rea- 
son. When the bulge Sersic index cannot be constrained 
due to spatial resolution and/or signal-to-noise limita- 

^ http://cas.sdss.org/dr7/en/tools/explore/obj.asp 



tions, its posterior probability distribution (which is fully 
mapped by GIM2D) will be uniformly flat between the 
minimum and maximum allowed values. Our allowed 
range of values was 0.5-8 based on previous studies of 
the Sersic index of spheroids. The median value of a flat 
posterior probability distribution (which we take to be 
the best-fit value) will therefore be around 4. Indeed, if 
we re-plot Figure [T5l for the entire sample with no i^-test 
selection, we see a very strong peak at Uf, = 4, but this 
peak reflects a lack of constraint on bulge profile rather 
than its actual shape. Third, there is an upturn in the 
Hh distribution at ni, > 7.5. We again examined a plot of 
Te versus e for these galaxies and found peaks at e ~ 
and e ~ 0.7 with a uniform distribution in sizes over the 
range r^ < 3" . Visual inspection of galaxies at e '-^ 
and e ^ 0.7 with smaller sizes {r^ ^ 0;'4) showed them 
to have nuclear, on-cenier sources. Galaxies at e ~ and 
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larger sizes (re ^ 2—3" ) did not exhibit any distinguish- 
ing characteristic as a group, and galaxies at e '^ 0.7 
and larger sizes (re ~ 2—3" ) had a bar+point source 
configuration. Galaxies with low and high n& values in 
Figure [T5] also have low and high Ug values in Figure [T31 
All of these different sub-classes of objects highlight the 
power of comparing different fitting models to identify 
different types of galaxy substructures. 

Galaxy half-light radii from nf, = 4 and free rib decom- 
positions {rfii^ni and r^iju respectively) are quite consis- 
tent over the full range of allowed n^ values (Figure fT6|). 
Not surprisingly, galaxy half-light radii differ between 
pure Sersic {rhi,ps) and bulge+disk fits at rig > 4 with 
rhi,pS being 50% larger at Ug = 8. This is entirely due 
to the fact that the half-light radii are calculated by in- 
tegrating best-fit models with different profiles in their 
outer wings. The bulge radii exhibit the expected shape 
from the well-known and strong covariance between rif, 
and re , and the choice of Uh will obviously have a signif- 
icant impact on the measurements of re- On the other 
hand, it is very important to note that the disk scale 
length does not appear to be affected by the Sersic index 
of the bulge for the majority (80%) of the galaxies in 
our sample, and that the scatter in the disk luminosity- 
size relation (Figure ^ therefore does not depend on the 
choice of bulge Sersic index. This apparent lack of de- 
pendence of rd on rib may be due to the lack of con- 
straint on bulge profile shape discussed earlier, but it 
may also be due to bulges being usually more compact 
than their disks. The covariance between measured bulge 
and disk parameters will be weaker in galaxies where the 
two components are more spatially distinct. More de- 
tails on the disk luminosity-size distribution are given in 
Simard 2011, in preparation. 

4.2.3. Sersic Model Companson with NYU Value-Added 

Catalog 

T he New York Universi ty Value-Added Galaxy Cata- 
log (" Blanton et al.l [2005al) provides Sersic model struc- 
tural parameters for galaxies in the SDSS spectroscopic 
sample. The details and tests of the NYU Sersic mea- 
surements including artifici al galaxy simulations a re de- 
scribed in the appendix of iBlanton et al.l ()2005b[) . We 
matched objects in our pure Sersic structural catalog 
with objects in the NYU catalog using MJD, PLATEID 
and FIBERID for cross-identifications, and a match was 
found for 666,740 objects. 

Figure [T7] shows the comparison between NYU and 
GIM2D galaxy Sersic half-light radii and indices. The 
trend in half-light radius shown in the left-hand panel 
is fully consistent with the NYU simulations if ng_gim2d 
and rhi,gim2d are taken to be equivalent to the input 
("true") values used for the NYU simulations (ni„ and 
''"so.ira respectively). The galaxy half-light radii from the 
NYU fits are smaller by about 20% than both their in- 
put simulation values and the GIM2D values for objects 
with rig > 5. There is an offset An ^ 0.3-0.4 between 
GIM2D and NYU Sersic indices at n^ = 1. This off- 
set depends on galaxy ellipticity: it increases from 0.2 at 
low e (< 0.1) to 0.5 at higher e (> 0.4). This dependence 
on ellipticity comes from the fact that the NYU profile 
fits were done on one-dimensional profiles extracted from 
two-dimensional images using circular annuli. The offset 
does not completely disappear even at low e possibly as a 



result of the fact that the NYU fits were done on r— band 
images only whereas our fits were done simultaneously on 
g— and r— band images, and a redder band will be more 
dominated by the redder, spheroidal (i.e., higher Sersic 
index) component of a galaxy. 

The comparison between the NYU and GIM2D Sersic 
parameters therefore shows good agreement given the 
differences in how the two sets of parameters were mea- 
sured. 

4.3. Data Tables and Some Cautionary Notes 



The data quality metrics used in Section [4 . 1 1 show that 
the GIM2D "SIM-|-SEXTDEBL-fGM2DBKG" dataset 
gives the most robust photometric results. The photo- 
metric data for this dataset are given in Tables [U [5] and [3] 
for the rif, = 4 bulge + disk, the free-n;, bulge -t- disk and 
the pure Sersic decompositions respectively. Two sets 
of galaxy half-light radii are given in these tables. The 
semi-major half-light radius Rhi of a galaxy was calcu- 
lated by individually collapsing the bulge and disk com- 
ponents onto their respective major axes, adding these 
two one-dimensional profiles into a global galaxy profile 
and computing the half light radius of this summed one- 
dimensional profile. The circular half light radius Rchi 
was computed by performing curve-of-growth photome- 
try in circular apertures on the intrinsic (i.e., not PSF 
convolved) GIM2D best-fit model image of the galaxy. 
The relationship between the two kinds of half-light ra- 
dius is easy to understand for a single component galaxy 
because it is given in this case by Rcm ~ RhiV^ — e where 
e is the ellipticity of this component. 

The best way to use the data listed in the structural 
parameter tables is to consider bulges and disks as sepa- 
rate galaxy sub-populations overlapping on the sky. As 
one would expect, bulge structural parameters are more 
reliable for brighter bulges (Figure ITSl) , and disk struc- 
tural parameters are more reliable for brighter disks (Fig- 
ure [T^ . The reliability of the structural parameters of a 
given subcomponent is largely independent of the other 
subcomponent (e.g.. Figure [T6l) . One should not study 
bulge or disk properties on the basis of a simple selection 
cut on bulge fraction. The bulge of a bright, low B/T 
galaxy can still be brighter than the bulge of a faint, 
high B/T galaxy. Bulge and disk subsamples should be 
selected on the basis of their magnitudes. 

Another important note of caution is related to the 
use of bulge fraction cuts to select early-type galaxy 
subsamples. As many previous studies have shown 
dim et al.l [200l iMcIntosh et al.l [2002t iTran et all [20031: 
iBlakeslee et'aLl boOGi: iSimard et al.l I2009D . earlv-tvpe 
galaxies should be selected using both bulge fraction and 
image smoothness. For example, a nuclear starburst in a 
relatively irregular galaxy would yield a high bulge frac- 
tion, and such a galaxy would be erroneously classified as 
an early- type if bulge fraction were the sole selection cri- 
terion. The agreement between visual and quantitative 
classification of galaxies has bee n shown to be exce llent 
when both parameters are used ([Simard et al.ll2009l ). 

Finally, internal dust should also be considered as a 
potential source of bias when selecting samples of bulges 
and/or disks for st udy. Disks are known to be dusty (e.g., 
[Driver et al.ll2007l and references therein), and this dust 
can broaden the intrinsic disk/bulge luminosity-size or 
luminosity-color distributions. One can deal with this 
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Figure 14. r-band GIM2D galaxy bulge fractions and Sersic indices from the tliree different fitting models. Top left: Distribution of 
the galaxy Sersic index ng from tlie single component, pure Sersic fits. Top right: Bulge fraction {B/T)„4 from the n;, = 4 bulge + disk 
decompositions versus galaxy Sersic index ng. Bottom left: Bulge fraction {B/T)f„ from the free ni, bulge + disk decompositions versus 
galaxy Sersic index ng. Bottom right: The difference between {B/T)f„ and {B/T)„4, as a function of bulge Sersic index nj. The greyscale 
represents the two-dimensional distribution of galaxies normalized by the total number of galaxies. The bin sizes for the two-dimensional 
distributions were A(ng) = 0.1 and A(B/T) = 0.02. The limits for the greyscale go from zero to 20% of the peak value of the normalized 
distribution. 



bias by either selecting galaxies with face-on disks or 
applying empirical internal dust extinction corrections. 
These corrections can be difficult to characterize, but 
the large number of galaxies included in the catalogs pre- 
sented here offers the opportunity of deriving these cor- 
rections as needed by various science programs (Simard 
2011, in preparation). 

5. SUMMARY 

We have performed bulge+disk decompositions for 
1.12 million galaxies in the Legacy area of the SDSS Data 
Release 7. Four decomposition procedures were used, 
and one of these procedures clearly produced more robust 
structural parameters, magnitude and colors when look- 
ing at three science-based data quality assurance metrics. 
The most reliable procedure included the following three 
important steps: 

• GIM2D-based sky background level determination 



• SExtractor object deblending 

• Simultaneous bulge+disk decomposition in g and r 

We also used three different fitting models: a. nt, — A 
bulge + disk model, a free-nt, bulge -f disk model and 
a pure Sersic model, and we provide a detailed compar- 
ison between the measured structural parameters from 
these fitting models with a mean to select the appropri- 
ate model for a given galaxy using the F-statistic. This 
comparison highlights the importance of this selection 
for a given science goal. Pure Sersic model fits might be 
better suited to studies of global galaxy colours whereas 
a cut on the i^-test probability Pps to select galaxies for 
which a bulge+disk model was required would be bet- 
ter for studies of bulge/disk colours. Using a low cut on 
both Pps and Pps probabilities would select relatively 
small but very robust samples of bulge+disk decompo- 
sitions. Full catalogs of structural parameters are in- 
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Figure 15. Distribution of GIM2D bulge Sersic index rij, from the 
free nj bulge + disk decompositions for galaxies with {B/T)f„ > 
and Pps < 0.32 and P„4 < 0.32. This distribution includes 52,897 
galaxies. 

eluded here with important cautionary notes on how to 
select bulge and/or disk subsamples, how to select differ- 
ent morphological types, the issue of bars and the issue 
of internal dust. These catalogs should provide an ex- 
tensive comparison set for a wide range of observational 
and theoretical studies of galaxies. 
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Figure 16. Comparison of r-band GIM2D radii measured from three different fitting models. Top left: Ratio of the galaxy half-light radii 
fhljn and r^j „4 from the free Jij, bulge + disk and ni, = 4 bulge + disk decompositions versus bulge Scrsic index ni,. Top right: Ratio 
of the galaxy half-light radii r^ipg and rf^i^i from the single component, pure Sersic fits and nj = 4 bulge -|- disk decompositions versus 
the galaxy Sersic index ng. Bottom left: Ratio of the bulge effective radii r^ f„ and r^ „4 from the free n;, bulge -|- disk and Uf, = 4 bulge 
-I- disk decompositions versus bulge Sersic index n;,. Bottom right: Ratio of the disk scale lengths r^f„ and r^j „4 from free n;, bulge -|- 
disk and Ji;, = 4 bulge -|- disk decompositions versus bulge Sersic index ni,. The greyscale represents the two-dimensional distribution of 
galaxies normalized by the total number of galaxies. The x-axis and y-axis bin sizes for the two-dimensional distributions were 0.1 and 
0.02 respectively. The limits for the greyscale go from zero to 50% of the peak value of the normalized distribution. 
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Figure 17. Comparison between NYU and GIM2D pure Sersic structural parameters. Left: Ratio of the NYU and circularized GIM2D 
galaxy half-light radii rso, sersic and rf^igi^2d versus the GIM2D galaxy Sersic index rig gi^2d- Right: Difference between the NYU and 
GIM2D galaxy Sersic indices versus GIM2D galaxy Sersic index ng^gim2d- The greyscale represents the two-dimensional distribution of 
galaxies normalized by the total number of galaxies. The x-axis and j;-axis bin sizes for the two-dimensional distributions were 0.1 and 
0.05 respectively. The limits for the greyscale go from zero to 30% of the peak value of the normalized distribution. 
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Figure 18. GIM2D bulge g and r photo metr ic e rror s as a function of bulge magnitude in n=4 bulge+disk decompositions. Bulge 
magnitudes were calculated using Equations I3cl and I3dl Errors on total galaxy magnitudes and bulge fractions were propagated through 
these equations to obtain the photometric errors. The data points are median values, and the error bars are the 16th and 84th percentile 
values. 
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Figure 19. GI M2D disk g and r photometric errors as a function of disk magnitude in n=4 bulge + disk. Disk magnitudes were calculated 
using Equations I3el and [M Errors on total galaxy magnitudes and bulge fractions were propagated through these equations to obtain the 
photometric errors. The data points are median values, and the error bars are the 16th and 84th percentile values. 
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SDSS structural parameters from ni,=4 bulge 



Table 1 

disk decompositions (Table available in online electronic version) 



Column 
Name 



Description 



ObjID 

z 

SpecClass 

Scale 

y-max 

9g2d 

rg2d 

9g2d,f 

rg2d,f 

A (fiber color) 

{B/T)g 

{B/T)r 

iB/T)gJ 

{B/T)rj 

Rhl,g 
Rhl.r 
Rchl,g 

Re 



<Pb 
Rd 



<Pd 

{dx)g 

{dy)g 

(dx)r 
(dy)r 

S2g 

Mg,g 
Mg,b 

Mr,g 
Mr,b 
Mr,d 

rib 
PpS 



SDSS Object ID 

SDSS Redshift (Spectroscopic if available. Photometric otherwise) 

SDSS SpecClass value (set to —1 if z is photometric or —2 if no redshift available at all) 

Physical scale in arcsec/kpc at redshift z 

Galaxy volume correction in Mpc^ (Equation [Jll 

g-band apparent magnitude of GIM2D output B+D model (Equation [ij 

r-band apparent magnitude of GIM2D output B+D model (Equation [TJ 

g-band apparent fiber magnitude of output B+D model 

r-band apparent fiber magnitude of output B+D model 

Delta fiber color defined as (g — T)gi^2d, fiber ~ id ~ r')sDSS, fiber (set to —99.99 if no SDSS fiber magnitudes available) 

g-band bulge fraction 

r-band bulge fraction 

g-band fiber bulge fraction 

r-band fiber bulge fraction 

g-band galaxy semi-major axis, half-light radius in kiloparsecs 

r-band galaxy semi-major axis, half-light radius in kiloparsecs 

g-hand galaxy circular half-light radius in kiloparsecs 

r-band galaxy circular half-light radius in kiloparsecs 

Bulge semi-major effective radius in kiloparcsecs (Equation [21l 

Bulge ellipticity (e = 1 — b/a, e = for a circular bulge) 

Bulge position angle in degrees (measured clockwrise from the +y axis of SDSS images) 

Exponential disk scale length in kiloparsecs (Equation [JJ 

Disk inclination angle in degrees (i = for a face-on disk) 

Disk position angle in degrees (measured clockwise from the +y axis of SDSS images) 

B+D model center offset from column position given by colc_g on SDSS corrected g-band image (arcsec) 

B+D model center offset from row position given by rowc_g on SDSS corrected g-band image (arcsec) 

B+D model center offset from column position given by colc_r on SDSS corrected r-band image (arcsec) 

B+D model center offset from row position given by rowc_r on SDSS corrected r-band image (arcsec) 

g-hand image smoothness parameter (as defined in Simard et al. (2009)) 

r-band image smoothness parameter (as defined in ^iniar^_et_al^ ^009)) 

g-hand GIM2D galaxy rest-frame, absolute magnitude (Equation l3all 

g-band GIM2D bulge rest-frame, absolute magnitude (Eauation l3cl l 

g-band GIM2D disk rest-frame, absolute magnitude (Equation I3c[l 

r-band GIM2D galaxy rest-frame, absolute magnitude (Equatio n I3bl l 

r-band GIM2D bulge rest-frame, absolute magnitude (Equation |3dll 

r-band GIM2D disk rest-frame, absolute magnitude (Equation [3ft 

Bulge Sersic index 

F-test probability that a B+D model is not required compared to a pure Sersic model 
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SDSS structural parameters from free n;, bulge 



Table 2 

disk decompositions (Table available in online electronic version) 



Column 
Name 



Description 



ObjID 

z 

SpecClass 

Scale 

y-max 

9g2d 

rg2d 

9g2d,f 

rg2d,f 

A (fiber color) 

{B/T)g 

{B/T)r 

iB/T)gJ 

{B/T)rj 

Rhl,g 
Rhl.r 
Rchl,g 

Re 



<Pb 
Rd 



<Pd 

{dx)g 

{dy)g 

(dx)r 
(dy)r 

S2g 

Mg,g 
Mg,b 

Mr,g 
Mr,b 
Mr,d 

rib 
PpS 

Pni 



SDSS Object ID 

SDSS Redshift (Spectroscopic if available. Photometric otherwise) 

SDSS SpecClass value (set to —1 if z is photometric or —2 if no redshift available at all) 

Physical scale in arcsec/kpc at redshift z 

Galaxy volume correction in Mpc^ (Equation [Jll 

g-band apparent magnitude of GIM2D output B+D model (Equation [ij 

r-band apparent magnitude of GIM2D output B+D model (Equation [TJ 

g-band apparent fiber magnitude of output B+D model 

r-band apparent fiber magnitude of output B+D model 

Delta fiber color defined as (g — T)gi^2d, fiber ~ id ~ r')sDSS, fiber (set to —99.99 if no SDSS fiber magnitudes available) 

g-band bulge fraction 

r-band bulge fraction 

g-band fiber bulge fraction 

r-band fiber bulge fraction 

g-band galaxy semi-major axis, half-light radius in kiloparsecs 

r-band galaxy semi-major axis, half-light radius in kiloparsecs 

g-hand galaxy circular half-light radius in kiloparsecs 

r-band galaxy circular half-light radius in kiloparsecs 

Bulge semi-major effective radius in kiloparcsecs (Equation [21l 

Bulge ellipticity (e = 1 — b/a, e = for a circular bulge) 

Bulge position angle in degrees (measured clockwrise from the +y axis of SDSS images) 

Exponential disk scale length in kiloparsecs (Equation [JJ 

Disk inclination angle in degrees (i = for a face-on disk) 

Disk position angle in degrees (measured clockwise from the +y axis of SDSS images) 

B+D model center offset from column position given by colc_g on SDSS corrected g-band image (arcsec) 

B+D model center offset from row position given by rowc_g on SDSS corrected g-band image (arcsec) 

B+D model center offset from column position given by colc_r on SDSS corrected r-band image (arcsec) 

B+D model center offset from row position given by rowc_r on SDSS corrected r-band image (arcsec) 

g-hand image smoothness parameter (as defined in Simard et al. (2009)) 

r-band image smoothness parameter (as defined in ^imar^_et_al^ ^009)) 

g-hand GIM2D galaxy rest-frame, absolute magnitude (Equation l3all 

g-band GIM2D bulge rest-frame, absolute magnitude (Eauation l3cl l 

g-band GIM2D disk rest-frame, absolute magnitude (Equation I3c[l 

r-band GIM2D galaxy rest-frame, absolute magnitude (Equatio n I3bl l 

r-band GIM2D bulge rest-frame, absolute magnitude (Equation |3dll 

r-band GIM2D disk rest-frame, absolute magnitude (Equation [3ft 

Bulge Sersic index 

F-test probability that a B+D model is not required compared to a pure Sersic model 

F-test probability that a free rib B+D model is not required compared to a fixed ni,=4 B+D model 
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Table 3 

SDSS structural parameters from pure Scrsic decompositions (Table available in online electronic version) 



Column 
Name 



Description 



ObjID 

z 

SpecClass 

Scale 

y-max 

9g2d 

rg2d 

9g2d,f 

^92d,f 

A (fiber color) 

Rhl,g 

Rhl.r 

Rchl,g 

Rchl.r 

e 

fpb 

{dx)g 

{dy)g 

(dx)r 

{dy)r 

S2g 

^9,9 
Mr,g 



SDSS Object ID 

SDSS Redshift (Spectroscopic if available. Photometric otherwise) 

SDSS SpecClass value (set to —1 if z is photometric or —2 if no redshift available at all)) 

Physical scale in arcsec/kpc at redshift z 

Galaxy volume correction in Mpc^ (Equation [Jll 

g-hand apparent magnitude of GIM2D output pure Sersic model (Equation [TJ 

r-band apparent magnitude of GIM2D output pure Sersic model (Equation [TJ 

g-band apparent fiber magnitude of output pure Sersic model 

r-band apparent fiber magnitude of output pure Sersic model 

Delta fiber color defined as (g — r)gi^2djiher ~ {g ~ ^)sDSS, fiber (set to —99.99 if no SDSS fiber magnitudes available) 

g-band galaxy semi-major axis, half-light radius in kiloparsecs 

r-band galaxy semi-major axis, half-light radius in kiloparsecs 

g-band galaxy circular half-light radius in kiloparsecs 

r-band galaxy circular half-light radius in kiloparsecs 

Galaxy ellipticity (e = 1 — b/a, e = for a circular galaxy) 

Galaxy position angle in degrees (measured clockwise from the +y axis of SDSS images) 

Pure Sersic model center offset from column position given by colc_g on SDSS corrected g-band image (arcsec) 

Pure Sersic model center offset from row position given by rowc_g on SDSS corrected g-band image (arcsec) 

Pure Sersic model center offset from column position given by colc_r on SDSS corrected r-band image (arcsec) 

Pure Sersic model center offset from row position given by rowc_r on SDSS corrected r-band image (arcsec) 

g-band image smoothness parameter (as defined in _Simar^_et_^l^ (2009)) 

r-band image smoothness parameter (as defined in Simard et al. (2009 )) 

g-band GIM2D galaxy rest-frame, absolute magnitude (Equation ISaf 

r-band GIM2D galaxy rest-frame, absolute magnitude (Equation I3bl l 

Galaxy Sersic index 



